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Abstract 



The cosmological relic density of the hghtest supersymmetric particle 
of the minimal supersymmetric standard model is calculated under the 
assumption of gauge and Yukawa coupling unification. We employ radia- 
tive electroweak breaking with universal boundary conditions from gravity- 
mediated supersymmetry breaking. Coannihilation of the lightest super- 
symmetric particle, which turns out to be an almost pure bino, with the 
next-to-lightest supersymmetric particle (the lightest stau) is crucial for 
reducing its relic density to an acceptable level. Agreement with the mixed 
or the pure cold (in the presence of a nonzero cosmological constant) dark 
matter scenarios for large scale structure formation in the universe requires 
that the lightest stau mass is about 2 — 8% larger than the bino mass, which 
can be as low as 222 GeV. The smallest allowed value of the lightest stau 
mass turns out to be about 232 GeV. 
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I. INTRODUCTION 



It is by now clear U| that, in a universe with zero cosmological constant, a combi- 
nation of cold plus hot dark matter is needed for fitting the data on cosmic microwave 
background (CMB) anisotropies and large scale structure in the universe, especially 
for essentially fiat spectrum of the primordial density fiuctuations. The energy density 
p of the universe is taken equal to its critical value pc (^2 = p/ pc = 1), as suggested by 
infiationary cosmology, and assumed to consist solely of matter [Vim = !)• About 10% 
of matter is baryonic {VLb ~ 0.1), while the rest (dark matter) contains a hot component 
with density equal to about 20% of the critical density {VLhdm ~ 0.2) and a cold one with 
^CDM ~ 0.7. The present value of the Hubble parameter in units of 100 km sec~^ Mpc~^ 
is taken to be /i ~ 0.5. Hot dark matter may consist of light neutrinos. This is compatible 
with the atmospheric and solar neutrino oscillations, within a three neutrino scheme, 
only if light neutrino masses are almost degenerate. A consistent super symmetric infia- 
tionary model with degenerate light neutrino masses providing the hot dark matter in 
the universe has been constructed in Ref . ||^ . Cold dark matter, in the case of vanishing 
cosmological constant, must satisfy the relation Qcdm h"^ ~ 0.175. 

Recent observational developments, however, seem to hint towards an alternative 
picture for the composition of the energy density of the universe with a nonvanishing 
contribution from something like a cosmological constant. Measurements of the cluster 
baryon fraction combined with the low deuterium abundance constraint on the baryon 
asymmetry of the universe, VLb h? ^ 0.02, suggest that the matter density is around 35% 
of the critical density of the universe {VLm ~ 0.35). Also, recent observations 0] favor the 
existence of a cosmological constant, whose contribution to the energy density can be as 
large as 65% of the critical density (I^a ^ 0.65), driving the total energy density close to 
its critical value as required by infiation. The assumption that dark matter contains only a 
cold component leads then to a 'good' fit of the CMB radiation and both the large scale 
structure and age of the universe data. Higher values of the Hubble constant {h ^ 0.65) 
are, however, required and, thus, ^cdm ~ 0.3. Moreover, the possibility of improving 
this fit by adding light neutrinos as hot dark matter appears to be rather limited. We 
can, thus, assume hierarchical neutrino masses in this case. A consistent supersymmetric 
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picture leading 'naturally' to hybrid inflation and employing hierarchical neutrino masses 
has been presented in Ref. |T^. In the presence of a nonvanishing cosmological constant, 
cold dark matter must satisfy flcDM h"^ ~ 0.125. 

Both these cosmological models with zero/nonzero cosmological constant, which pro- 
vide the best flts to all the available data, are equally plausible alternatives for the 
composition of the energy density of the universe. Thus, taking into account the obser- 
vational uncertainties, we will restrict ^cdm h"^ in the range 0.09 — 0.22. 

The lightest supersymmetric particle (LSP) of the minimal supersymmetric standard 



model (MSSM) is one of the most promising candidates for cold dark matter ||TT],|T2[. This 
is normally the lightest neutralino and its stability is guaranteed by the presence of a dis- 
crete Z2 matter parity, which implies that supersjTiimetric particles can disappear only 
by annihilating in pairs. The cosmological relic density of the lightest neutralino can be 
reliably computed, for various values of the parameters of MSSM, under the assumptions 
of gauge coupling uniflcation and radiative electroweak breaking with universal bound- 
ary conditions from gravity-mediated supersymmetry breaking (see e.g., Refs. [|T^-[T5|). 
Coannihilation |]T6[ of the LSP with the next-to-lightest supersymmetric particle (NLSP) 
turns out to be crucial in many cases [0,^,0 . 

The purpose of this paper is to estimate the lightest neutralino relic density in a 
speciflc MSSM framework [|l8l of the above variety, where the three Yukawa couplings of 
the third family of quarks and leptons unify 'asymptotically' (i.e., at the grand unifled 
theory (GUT) scale Mqut ~ 10^^ GeV). This can arise by embedding MSSM in a 
supersymmetric GUT based on a gauge group such as 5*0(10) or Eq, where all the 
particles of one family belong to a single representation. It is then obvious that requiring 
the masses of the third family fermions to arise primarily from their unique Yukawa 
coupling to a particular superfleld representation (say a 10-plet of 5*0(10)) predominantly 
containing the electroweak higgses guarantees the desired Yukawa coupling uniflcation. 
This scheme predicts large tan/5 mt/mi), as well as the successful 'asymptotic' mass 
relation = rrih. The supersymmetric particle spectrum, top quark mass and higgs 
scalar masses in this model have been studied in Refs. p9|-pT|. The top quark mass is 
'naturally' restricted to large values compatible with the present experimental data and 
the supersymmetric particle masses are predicted relatively large. The lightest neutralino 
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is an almost pure bino, whereas the NLSP is the hghtest stau mass eigenstate. 

Coannihilation of the bino with the NLSP turns out to be of crucial importance for 
keeping the bino relic density at an acceptably low level. This implies that the lightest 
stau must not be much heavier than the bino so that coannihilation can be effective. 
Moreover, increasing the lightest stau to bino mass ratio leads to a larger bino mass which 
further enhances its relic density. Lightest stau masses of about 2 — 8% larger than the 
bino mass are required for obtaining Qcdm h"^ in the range 0.09 — 0.22. It is interesting to 
note that, for smaller 'relative' mass gaps between the lightest stau and the bino, ^cdm h'^ 
rapidly decreases and becomes unacceptably small. The values of this mass gap which we 
find here combined with the fact that the bino mass turns out to be greater than about 
222 GeV make the lightest stau a phenomenologically interesting charged sparticle with 
mass which can be as low as ~ 232 GeV. Our analysis provides quite strong restrictions 
on the sparticle spectrum of MSSM with Yukawa coupling unification. 

In Sec.0, the MSSM with Yukawa coupling unification is introduced and its parame- 
ters and sparticle spectrum are constrained. In Sec.|T|, the relic LSP (lightest neutralino) 
density is calculated by taking into account its coannihilation with the NLSP (lightest 
stau). In particular, the bino annihilation cross section is estimated in Sec. Ill A] , whereas 
Sec. [lllB] is devoted to the evaluation of the relevant coannihilation cross sections. Our re- 
sults on Q LSP h"^ are presented and their consequences are discussed in Sec. [III Q . Finally, 



our conclusions are summarized in Sec. IV. 



II. MSSM WITH YUKAWA UNIFICATION 



We consider the MSSM embedded in some general supersymmetric GUT based on a 
gauge group such as 5*0(10) or (where all the particles of one family belong to a single 
representation) with the additional requirement that the top, bottom and tau Yukawa 
couphngs unify at the GUT scale Mqut- This requirement is easily guaranteed by 
ensuring that the masses of the third family fermions arise primarily from their unique 
Yukawa coupling to a single superfield representation which predominantly contains the 
electroweak higgses. We further assume that the GUT gauge symmetry breaking occurs 
in one step. Ignoring the Yukawa couplings of the first and second generation, the effective 
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superpotential below Mqut is 

W = eiji-htW^Qit' + h,HlQib' + hrHlLir' + fiHlHi) , (1) 

where = {t, b) and L3 = (z/^, r) arc the quark and lepton SU{2)l doublet left handed 
superfields of the third generation and t^, If and r'^ the corresponding SU{2)i singlets. 
Also, i^i, H2 are the electroweak higgs superfields and €12 = +1. The gravity-mediated 
soft supersymmetry breaking terms in the scalar potential are given by 

a,b 

[eiji-AthtH'^Qif + AbhbHiQib'' + ArKHiLir + BptHlHi) + h.c.) , (2) 

where the 's are the (complex) scalar fields and tildes denote superpartners. The 
gaugino mass terms in the Lagrangian are 

1{M,BB + M2 ^ WrWr + Ms ^ ~ga~ga + h.c.) , (3) 

^ r=l a=l 

where B, Wr and cja are the bino, winos and gluinos respectively. 'Asymptotic' Yukawa 
coupling unification implies 

htiMcuT) = K{Mgut) = K{Mgut) = ho . (4) 

Based on = 1 supergravity, we take universal soft supersymmetry breaking terms at 
Mqut, i-e., a common mass for the scalar fields mo, a common trilinear scalar coupling 
Aq and Bq = Aq — mo. Also, a common gaugino mass M1/2 is assumed at Mqut- 
Our effective theory below Mqut depends on the parameters (/^o = 1^{Mgut)) 

mo, Ml/2, Ao, iio, aa, Mqut, ho, tan/3 . 

The quantities ao — 9%!^'^ (Qg being the GUT gauge coupling constant) and Mqut are 
evaluated consistently with the experimental values of ctem, cts and siv? 6w at mz- We 
integrate numerically the renormahzation group equations (RGEs) for the MSSM at two 
loops in the gauge and Yukawa couplings from Mqut down to a common supersymmetry 
threshold Ms ~ 1 TeV. From this energy to mz, the RGEs of the nonsupersymmetric 
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standard model are used. The set of RGEs needed for our computation can be found in 
many references (see, for example, Ref. |^). We take ^^(m^) = 0.12 ± 0.001 which, as 
it turns out, leads to gauge coupling unification at Mqut with an accuracy better than 
0.1%. This allows us to assume an exact unification once the appropriate supersymmetric 
particle thresholds are taken into account. Our integration procedure relies on iterative 
runs of the RGEs from Mqut to low energies and back, for every set of values of the 
input parameters, until agreement with the experimental data is achieved. The value of 
tan/? at Ms is estimated using the experimental input mr{mr) = 1.777 GeV and Ms is 
fixed to be 1 TeV throughout our calculation. Assuming radiative electroweak symmetry 
breaking, we can express the values of the parameters /i (up to its sign) and B at Ms in 
terms of the other input parameters by means of the appropriate conditions 



_ m^^ - m^^ tan P 1 _ 2Bf^ 

^ ~ tan2 p-l ^ra,,smzp- ^ ^ , (5j 



where m^j, are the soft supersymmetry breaking scalar higgs masses. Here, fol- 



lowing Ref. [23|, we used the tree-level renormalization group improved scalar potential 
minimized at a scale comparable to the mass of the stop quark. This is adequate for our 
purposes since, as we find, the corrections to /i from the full one-loop effective potential 



in Ref. [2^] are negligible. The sign of /i is taken to be negative (with our conventions), 
which leads to acceptable predictions for 6 — 57 in models with large tan/5 [0. 

The common value of the third generation Yukawa coupling at Mqut is found by 
fixing the top quark mass at the center of its experimental range, mt{'mt) = 166 GeV. 
The value obtained for mh{mz) after including supersymmetric corrections is somewhat 
higher than the experimental limit We are left with mo, M1/2 and Aq as free input 
parameters. Our results, as it turns out, depend very little on the exact value of Aq 
which is, thus, fixed to zero in our calculation. The values of mo and M1/2 are found as 
functions of the tree-level mass niA of the CP-odd higgs scalar A, for each 'relative' mass 
splitting between the NLSP (lightest stau) and the LSP (almost a pure bino), as we will 



explain later. The value of is evaluate at Ms which is comparable with y/fn^m^ 



Although the full one-loop corrections to (from Ref. |2^) are not totally negligible, 
we will ignore them here since their effect on the LSP relic density is small. 

The LSP is the lightest neutralino x- The mass matrix for the four neutralinos is 



5 



/ Ml 


-mzsw cos P 
\ mzSw sin (3 





M2 
mzcw cos /? 
—mzCw sin /? 



— m^svy cos j3 
mzCw cos /? 




m^SH/sin/? \ 
— mzCvt/ sin (3 







(6) 



in the {—iB, —iW^, Hi, H2) basis. Here sw = sin6'iy, cw = cos6'iy, and Mi, M2 are the 
mass parameters of B, W3 in Eq.(^). For the values of fi obtained from the radiative 
electroweak breaking conditions here (/i/Mi/2 ~ 1-2), the hghtest neutrahno turns out 
to be a bino, B, with purity > 98%. 

Large b and r Yukawa couphngs cause soft supersymmetry breaking masses of the 
third generation squarks and sleptons to run (at low energies) to lower physical values 
than the corresponding masses of the first and second generation. Furthermore, the large 
values of tan /? implied by the unification of the third generation Yukawa couplings lead 
to large off-diagonal mixings in the sbottom and stau mass-squared matrices. These 
effects make the physical mass of the lightest stau significantly lower than the masses of 
the other squarks and sleptons (see below). The NLSP is, thus, the lightest stau mass 
eigenstate f2 and its mass is obtained by diagonalizing the stau mass-squared matrix 
ml + m?^ + m|(— 1/2 + s^) cos 2/? m^(A^ + /xtan/?) 

m^(yl^ + /itan/3) + m?^ — m|s^ cos 2/? 

in the gauge basis (f^, tr). Here, rrir^^^^ is the soft supersymmetry breaking mass of 
fL(R) and uIt- the tau lepton mass. The stau mass eigenstates are 

T2/ \-S9 Cg) \fR 

where sg = sin 6', cg = cos 6, with 6 being the — mixing angle. Another effect of the 
large values of the h and r Yukawa couplings is the reduction of the mass of the CP-odd 
higgs boson and, consequently, the other higgs boson masses to smaller values. 

The authors of Ref. found that, for every value of and a fixed value of (m^), 
there is a pair of minimal values of tuq and M1/2 where the masses of the LSP and f2 are 
equal. This is understood from the dependence of on mg and M1/2 given in Ref. 
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m , 



aM\i2 — Pml — const. , 



(9) 



where all the coefficients are positive and a and /3, which depend only on mt{mt), are 
~ 0.1 (the constant turns out to be numerically close to Equating the masses of the 
LSP and f2 is equivalent to relating mg and Mi/2- Then, for every m^, a pair of values 
of mo and M1/2 is determined. Note that Eq.(|^) implies the existence of an upper bound 
on TTiA since m\ < aMf^^ ■ We set here an upper limit on M1/2 equal to 800 GeV, which 
keeps the sparticle masses below about 2 TeV consistently with our choice for M5 (=1 
TeV). This limit constrains m^i to be smaller than ^ 220 GeV. On the other hand, the 
experimental searches for the lightest CP-even neutral higgs boson h with mass rrih set 
a lower limit on ttla- Taking into account radiative corrections |]28|p9| in calculating rrih, 
we found that this lower limit on is about 95 GeV. The highest values of rrih, which 
are obtained as rriA increases to its upper limit, lie between 125 and 130 GeV. 



Following the procedure of Ref. [^, one can determine mo and M1/2 not only for 
equal masses of the LSP and NLSP but for any relation between these masses. We fix 
Tnt{mt) = 166 GeV (tan/3 ~ 52.9). For every m^ and a given 'relative' mass splitting 
= (mfj — my,)/my, between the NLSP and LSP, we find mo and M1/2. They are 
depicted in Fig.|l| as functions of m^ for Af2=0.02 and 0.08 (see Sec. [lll C| ). We observe 
that, for fixed tua, M1/2 increases with A^^. Thus, mo and the sparticle masses increase 
too with (see Eq.(^)). Also, for fixed M1/2, rnA is a decreasing function of Af^. 
As a consequence, the upper bound on niA (corresponding to M1/2 = 800 GeV) gets 
reduced as Af^ increases. This is why the curves in Fig.[l| which correspond to higher 
's terminate at smaller rriA 's. As we will see, the cosmological bounds on Qlsp h'^ 
will constrain Af^. The relevant part of the sparticle spectrum as a function of tjia, for 
Af2=0.047, is shown in Fig.^. The LSP mass, for Af2=0.02, is also included. 



III. LSP RELIC DENSITY 

We now turn to the calculation of the cosmological relic density of the lightest neu- 
tralino x (almost pure B) in MSSM with Yukawa unification. As mentioned in SecJI], 

h"^ increases to unacceptably high values as m^ becomes larger. Low values of rriy. are 
obtained when the NLSP (f2) is almost degenerate with x- Under these circumstances, 
coannihilation of x with f2 and is of crucial importance reducing further the x T^^^ic 
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density by a significant amount. The important role of coanniliilation of the LSP with 
sparticles carrying masses close to its mass in the calculation of the LSP relic density 
has been pointed out by many authors (see e.g., Refs. |T^,|T^JT^). Here, we will use the 
method described by Griest and Seckel [1^. Note that our analysis can be readily applied 
to any MSSM scheme where the LSP and NLSP are the bino and stau respectively. 
The relevant quantity, in our case, is the total number density 

n = n^ + nf^ + Uf* , (10) 

since the f2 's and f2 's decay into x 's after freeze-out. At cosmic temperatures relevant 
for freeze-out, the scattering rates of these (nonrelativistic) sparticles off particles in the 
thermal bath are much faster than their annihilation rates since the (relativistic) particles 
in the bath are considerably more abundant. Consequently, the number densities 
{i = X, ^2) proportional to their equilibrium values ra^^ to a good approximation, 
i.e., Ui/n ~ rf^ jn^'^ = ri. The Boltzmann equation (see e.g., Ref. pO]) is then written as 



- = -3Hn-{a^ffv){n'~{n^''r) , (11) 

where H is the Hubble parameter, v is the 'relative velocity' of the annihilating particles, 
(■ ■ ■) denotes thermal averaging and cxe// is the effective cross section defined by 

with (Tjj being the total cross section for particle i to annihilate with particle j averaged 
over initial spin and particle-antiparticle states. In our case, (Te// takes the form 

For Tj, we use the nonrelativistic approximation 

(?.(l + A.)3/^e-^-- 
ri{x) = , (14) 

9eff 

9eff{x) = Y.9^{l + A,f^e-^'\ A, = (m, - m^)/m^ . (15) 

i 

Here gi = 2, 1, 1 (i = x, f2, ) is the number of degrees of freedom of the particle species 
i with mass rrii and x = rrij^/T with T being the photon temperature. 
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In Table I, we list all the Feynman graphs included in the calculation of the effective 
cross section. The exchanged particles are indicated for each relevant pair of initial and 
final states. The symbols s{x), t{x) and u{x) denote tree-graphs in which the particle x 
is exchanged in the s-, t- or u-channel. The symbol c stands for 'contact' diagrams with 
all four external legs meeting at a vertex. H and denote the heaviest neutral and 
the charged higgs bosons, while e, en, u and d represent the first and second generation 
charged leptons, charged right handed sleptons, up- and down-type quarks. The other 
possible reactions f — > h[H]A, h[H]'y, h[H]Z, AZ, iJ^H^^ or i/z/ (i/ stands for all three 
neutrinos) have not been included since they are utterly suppressed by small couplings 
and/or heavy masses. Also, the tiny contributions from graphs with h and H exchange 
in the s-channel, in the cases of uu, dd, ee final states, are left out. Some of the graphs 



listed here have not been considered in previous works |14| with small tan/3. 



TABLE I. Feynman Diagrams 



Initial State 


Final State 


Diagrams 


XX 


TT 


t(fi,2), U{f 1^2) 




ee 


t{eR), u{eR) 


XT-2 


Th, tH, tZ 


sir), t(fi,2) 




tA 


sir), tin) 




TJ 


sir), t(f2) 




TT 


tix), uix) 




hh, hH, HH, ZZ 


sih), siH), t(fi,2), M(fi,2), c 




AA 


sih), siH), tin), uin), c 




H+H~, W+W 


sih), siH), 5(7), siZ), c 




77, iZ 


tir2), uif2), c 




tt, bb 


sih), siH), s(7), siZ) 




TT 


sih), siH), s(7), siZ), tix) 




UU, dd, ee 


s(7), siZ) 
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The relic abundance of the LSP at the present cosmic time can be calculated from 
the equation |TB|,|5D[ 

1.07 X 10^ GeV~i 



gl^'^Mp Xp" CTeff 



with 



eff = Xf {creffv)x ^dx . 



(16) 



(17) 



Xp 



Here Mp = 1.22 x 10 GeV is the Planck scale, ~ 81 is the effective number of 
massless degrees of freedom at freeze-out and xp = m^/Tp, with Tp being the 
freeze-out photon temperature calculated by solving iteratively the equation [p0| , plj 

0.038 Qeffixp) Mp (c + 2)c {aeffv){xp) 



Xp = In 



1/2 1/2 

9* ^p 



The constant c is chosen to be equal to 1/2 ||3T|. The freeze-out temperatures which we 
obtain here are of the order of m^jl'b and, thus, our nonrelativistic approximation (see 
Eq.(p^) is justified. Under these circumstances, the quantities OijV are well approximated 
by their Taylor expansion up to second order in the 'relative velocity'. 



The thermally averaged cross sections are then easily calculated 

3,3/2 



{aijv){x) 



2v^Jo 

Using Eqs.dll), (0), (^) and (ID, one obtains 



dvv'^{aijv)e ^'^ = aij + Qhij/x . 



(19) 



(20) 



where we sum over (ij) = (xx), (xra) and {f2f^*^) with a (,) = a^^f^ + a^^^*, 6 (.) 
6f2f2 + bf^f* and a(ij), given by 

f"^ dx 



(21) 



I Xp X 



ri{x)rj{x) , %) = Qc{^ij)Xp 



°° dx 



Xp X 



ri{x)rj{x) 



(22) 



Here C(,j) = 1, 4, 2 for (tj) = (xx), (^^2) and (^2^2^*^)- For A^^ = 0, = 1/4, 1/2, 1/8 
((u) = ixx), {Xf2), {r2ft^)), while %) = 3a^ij)/xp. 
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A. Annihilation cross section 



The fact that the LSP (x) is an almost pure B imphes that the main contribution 
to its annihilation cross section arises from sfermion (squark, slepton) exchange in the t- 
and u-channel leading to // final states (/ is a quark or lepton). The s-channel diagrams 
are suppressed since the values of obtained here are always far from 'mz/2 and mij2 
(see e.g., Ref. ||13[)- Moreover, diagrams with quarks in the final state are suppressed 
relative to the ones with leptons because of the heavier masses of the exchanged squarks 
and the smaller quark hypercharges. As mentioned in Sec.|I|, under the assumption of 
unification of the third family Yukawa couplings, mf^ is smaller than the masses of the 
other sleptons, hence the production of rf is enhanced relative to the production of 
lighter leptons. 



Using the partial wave expansion of Ref. [13] and neglecting the masses of the final 



state leptons, we evaluate the coefficients a^^ and h^^ in Eg . (p^ . They are found to be 



^ 2 2-1^2-1^2 ( ^ 



Si 



(23) 



hx 



^2 



+ 2 



m?(m^ + m|^) 



127rcf^ 



(24) 



m'^+m'^^ 2 and Eg 



where Y^n) = -l/2(-l) is the hypercharge of tl(r), 
with mg^ being the common (see below) mass of the right handed sleptons cr, jiR of the 
two lighter families. Some comments are now in order: 



i. The presence of a nonvanishing coefficient a^^^ is due to the large values of tan/3 
which lead to an enhancement of the off-diagonal terms in the stau mass-squared 
matrix in Eq.(|^). Indeed, this coefficient is negligible in the case of small — tr 
mixing (i.e., for low tan/5) where the f2 essentially coincides with f^. This is due 
to the fact that the s-wave contribution, which is the only contribution to a^^ , is 
suppressed by factors of the final state fermion mass as one can show by employing 
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Fermi statistics arguments []TT|. For large tan/3, however, this suppression is not 
complete and a^^ is proportional to sin^ 6. Despite the fact that a^^ is smaller than 
h^^ , its contribution to o"e// in Eq.(^T]) is of the same order of magnitude as the 
one of h^y. which enters in this equation divided by a relative factor < x^/S ~ 8 — 9. 

ii. The main contribution to h^^ arises from the first term in the bracket in the right 
hand side of Eq. (p^ ) . The second term in this bracket is due to — f/j mixing. 



iii. The last term in the right hand side of Eq.(|^) represents the contribution of the 
two lighter generations. Their right handed sleptons are considered degenerate with 
mass me^. The off-diagonal elements in the slepton mass-squared matrices of the 
lighter families are negligible. The values of mg^ are bigger than irtf^ and hence the 
corresponding contributions to hy^y^ are smaller than the ones from the f2 exchange. 
This is a major difference from models with low tan /5, where the contributions of 
all three diagrams with exchange of right handed sleptons are similar. 

iv. The contribution to h^^ of the diagram with a fi exchange is small and, although 
taken into account in the computation, is not displayed in Eq. (p^ . We find that this 
contribution is suppressed by about 1/6 — 1/8 compared to the contribution of each 
of the lightest generations. This can be understood by the following observation. 
Despite the fact that the values of the mass in the propagator of this diagram, irif^ , 
are not much higher than mg^, its main contribution contains a factor CqY^. 



B. Coannihilation Cross Sections 



The contributions of the various coannihilation processes listed in Table I to the 
coefficients aij and hij {ij ^ xx) in Eg . ([T9|) are calculated using techniques similar to the 
ones in Ref . . Leptons and quarks (except the t-quark) in final states or propagators 



are taken to be massless. On the contrary, the h and r Yukawa couplings are not ignored 
since, in our case where tan/3 is large, their influence turns out to be very significant. 
The most important contributions to o"e// in Eq.(pT]) arise from the aij 's in the case 
of coannihilation. In Table II, we list some of the processes contributing to the aij 's 
{ij 7^ xx) together with the analytical expressions for their contributions. 
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TABLE II. Contributions to the Coefficients Uij {ij ^ xx) 



Process 


Contribution to tfie Coefficient 


XJ2 - 




6^(1 - ■mlf{2YLYiighrT[^S0C0gh/{mf^ - rfilm^) - cos2%i/ 
(m?^ + m^{m^, - ml))] + [gl^^ + gl/{m^ml - mr,)\slYl + c^Fj) 

^-gUcffl + slYl)/{ml + m^(m,, - m^)^ - 25,c,(yi - Fj) 
9h9hi/{rhf2 - fh^rhl){m\ + m^(m^2 - m|))}/327rc^mf2(mf2 + m^) 


XJ2 - 


— >• TJ 


e-\slYl + c^Fj)/167rc^mf2(mx + m^J 


XJ2 - 


^tZ 


e2(l - m|){m^(l - ffiD^igl^zi^yl + c^^J)/(m|% - m^,)^ 

+9UMyi + + ^^x(^^2 - ^D? 

-'^9fif2z9f2f2zSeCe{Yl - Yl)/{mf^ - m%m^){m\ + m^{fhf^ - m|))] 
-2gz{fhl - ^T[9f27^z{LrslYl + R^clYi)/{fh^ - fhlfh^) 
-gf^f^zSeCe{LrYl - RrY^)/ {m\ + m^{rfif.^ - m|))] 
WziLlslYl + i??c2rj)(l + ml - 2m|)(l + m,)}/327rc^m| 


T2T2 


— > TT 


e\stYt + c^ri)m|/7rc^^E^ 
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^^Z 


-e^fi'^f2z(^l - 4)/167rm^ 




zz 


(1 - - 4) + l2g,gnzz9Uzmlml)/{rhl - 4) 

-4(??.(?/.zz^?|,f2z"^i(A - rnlPi)/{l + m?^ - ml){ml - A) + {h ^ H)] 
+^4f2z"^l^3/(m| - 2)2 + 2ghghzz9H9HzzPi/{ml - 4:){mjj - 4) 
-84f.z^Uz<[n - 3m?^(m| - 2)]/(l + - m|)(m| - 2) 
+4m|^<7|,.,z[< A + (1 - nil)'P2 - '2rhlP4/{l + m?^ - ml)'}/ 

647rm|m?2 


T2T2 


w+w- 


(1 - m^)^/^(4 - 4m^ + 3m^)[^hy;,H^+iy-/(m^ - 4) 

+9H9HW+w-/i^H - 4) + 5f^2^2H^+ty-m|JV327rm^m|2 


T2TI 


ti 


3(1 - rh'tnh9hu/{ml - 4) + ^^^H«/(m| - 4)]V47rm|^ 
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Here hat (or bar) over a quantity indicates that this quantity is measured in units of m^^ 
(or + nif^) and the g 's will be defined shortly. Also, L^- = I — 2s^, Rr = — and 



Pi(2)=3m|-(+)4m| + 4, P3 



3m^ 



iTTlr 



P4 = 3m| - 3m| - 4m| + 4 , P5 = 3m| - 5m| + 2 . 



(25) 



The contribution of the process x^2 tH (or tA) to the coefficient a^r2 is obtained 
from the expression for rh in Table II by replacing h hj H (or A and cos 26' by 

1). For the contribution to af^r* of each of the five processes with two higgses in the final 
state (see Table I), a general formula can be given: 

.1. 1 



A. 



4 -ml 



+ 



A 



H 



4 — mjj 



+ 



(4 



\2U/2 



(4 



v2\l/2 



4Ai 



+ 



4A, 



^ff-, + ^ff. — 2m|, — 2 ' mf^ + — 4 



Ac"^f2 , 



(26) 



where the Hp, Hg stand for h, H, A, H~^, H~ , the factor 1/2 enters only for identical 
particles in the final state and Xh, Xh, Ai, A2, Ac correspond to the diagrams s{h), s{H), 
^(^1,2) (or u{fi^2)), c in Table I and are shown in the Table III. 

Table III. The A Symbols 



Process 


Xh 


Xh 


Ai 


A2 


Ac 


^2^2 ~^ 


ghQhhh 


9H9hhH 


9I1 


9I 


9f2T2hh 


^2X2* hH 


gh9hhH 


9H9hHH 


9hl9Hl 


9h9H 


9f2f2hH 


f2T2 ^ HH 


ghghHH 


gHguHH 


9m 


9h 


9f2T2HH 


f2f2 A A 


ghghAA 


guguAA 


-9ai 





9f2T2AA 


f2f; H+H- 


9h9hH+H- 


9h9hh+h- 








9f2f2H+H' 



The g 's in Tables II and III correspond to vertices with the particles entering indicated 
as subscripts. The simplest ones are (for Feynman rules, see e.g., Ref. with 

gfif2Z = gz(-sece) , 5'f2?2^ = 9zisl - 2s^) , gf2f2W+w- = 9'^sl/2 , (27) 

where gz = g/2cw with g being the SU{2)l gauge coupling constant. Note that gA = 
9f2f2A = 0. The more complicated g 's are arranged in the Table IV. 
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Table IV. g Symbols 



g bymbol 


bxpression 


gh[H]i [ = gfif2h[H]) 


gzmz sm[- cosj [a + p) [Lr + RrjSece 
-\-gmr cos 26 (At sin[— cos] a — cosfsinja) /2mw cos (3 


9h[H] (= 9f2f2h[H]) 


—gzfnz sin[— cos] (a + P){LrSg — RtC^) — {gmr/mw cos (3) 
{—mr sin[— cosjo; — soCe{AT sin[— cosjct — //cos [sin] a)} 


„ I „ \ 

9a\ i= gfxf^A) 


gmr (At tan P — ij) / 2mw 


9f2f2hh[HH] 


— [+]9z cos2«(L^s^ — Rrcj) — g'^(sin[cos]a/ cos/5)^m^/2m^ 


9f2f2hH 


2 • o / r /o 2 1 2 /o 2 2 a\ 2 /n 

g sin2Q;(— Li-/2c^ + m1;./2my^ cos p)s^/2 
+5' sm zq;(— tan c^i^ + m^/2m^^ cos p)Cq/2 


9f2f2AA 


— gf^ COS 2p[LrS0 — KtCq) — 5^ tan p(mT-/myyj /z 


9f2f2H+H- 


g'^cos2/?((l - Lr/2c^r)sg - tan^ cyw'C|j/2 

„2+„„2/3/™ /™ \2„2 /o 

— tan pym-r/mw) Cg/Z 


9hhh[HHH] 


—3gzTnz sin[cos](a + P) cos 2a 


9hHH[hhH] 


fl'z^z{sin[cos] (a + /?) cos 2q; + 2 cos[— sin] (a + /?) sin 2q;} 


9h[H]AA 


—gzT^z sin[— cos] (a + /?) cos 2/3 


9h[H]H+H- 


—g{mw sin[cos] (j3 — a) + mz sin[— cos] (a + j3) cos 2j3 /2cw} 


9h[H]ZZ 


guiz sin [cos] (/3 — cc)/ 


9h[H]W+W- 


g^rnvK sin [cos] (/5 — a) 


9h[H\tt 


— gi(cos[sin]Q;/ sin/3)(mt/2miy) 


yhJJ\rT 


/y(sin[— eosjn / cos :7)(/;/^/2/;/ir) 


9Att 


—g tan P (mr/2mw) 



Here we have defined 

tan2Q; = tan2/3 (m^ + m|)/(m^-m|) , - 7r/2 < a < . (28) 

We do not show explicitly the small contributions to af^f* of the processes with 66 
and rf in the final state. They are, however, taken into account in the computation. 
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The contributions to a^;.^^* of the processes with uu, dd and ee in the final state vanish 
(these processes contribute only to h 's). Also, the coefficients hij {ij ^ xx)i although 
included in the calculation, are not displayed since their contribution to aeff is, in general, 
negligible. Note that many of the couplings and terms listed above have not been included 



in previous calculations |14] with small tan/3. Some comments are now in order: 



i. All five processes for the coannihilation of x with f2 listed in Table I give more 
or less comparable contributions to the coefficient a^f^ (the leading contribution 
comes, in general, from xf2 Th). The relative contribution of fe^fj to ^{xf2) in 
Eg. (|2T|) turns out to be essentially independent of (95 GeV< < 220 GeV). 
This contribution varies from about 5% to about 8% as increases from to 0.1 
(this is the relevant range of Afj as we shall see). 

ii. The major contributions to a-^-{-) come from the processes f2T2 hh, ti and 
^2T^2 TT. However, many of the other relevant processes in Table I (like 
^2T^2 ^ 77' HH, AA, H^H~, 7Z) have, in general, important contributions 
which cannot be neglected (T2f| ZZ, for large Af^ 's and 's, gives a major 
contribution). Also, the reaction f2f"2 —>■ hH {W~^W~) is enhanced for low values 
of rriA (and A~). The relative contribution of b. _(,) to a,, .m., which can be either 
positive or negative, is less than about 1% for all relevant values of parameters. 

iii. For = 0, the contribution of the x annihilation to o"e// is very small (~ 0.4%). 
The corresponding contributions of ^"(^fj) span the ranges 27 — 24% and 
73 — 76% respectively as rriA varies from 95 to 220 GeV. For Af^ = O-l; however, 
the annihilation of x 's becomes very significant accounting for about 33 — 31% of 
(Te//. The most important contribution 58% of CTeff), in this case, comes from 
the coannihilation of x with f2, whereas f2 coannihilation with f2 or accounts 
for about 9 — 11% of o"e//. We see that, although x annihilation is negligible for 
small Af2 's, it is strongly enhanced at higher values of Af^ . This is due to the fact 
that the abundance of f2 's decreases relative to the one of x 's as Af^ increases. 

iv. For Af2 = O5 the contributions of b^f^ and b^^^(,) to o"e// cancel each other partially 
and, thus, an accurate result (error ^ 0.5%) can be obtained by ignoring these b 's. 
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For =0.1, however, the contribution of fo^^f^ dominates strongly over the one 
of &-^-(*) and gives ^ 4 — 5% of (Xg//. Consequently, our results can be reproduced 
with an accuracy better than ^ 5% by using, for coannihilation, just the a 's. Their 
analytical expressions have been given earlier in this Section. On the contrary, the 
cannot be ignored since its contribution to o"(^^) can be as high as 80% and the 
annihilation of x 's is very significant at higher 's. 

C. Results on Qlsp 

We can now proceed to the evaluation of ^lsp h"^- The top quark mass mt{mt) is again 
fixed at 166 GeV. For given values of and m^, all the particle physics parameters of 
the model are determined (see Sec.|I|). The freeze-out parameter xp can then be found 
by solving Eq.(|T8|) and aeff is evaluated from Eq.(|21|). The LSP relic abundance VLy. is 
estimated using Eq.(|T6|) and is depicted in Fig.|^ as function of for Af2=0, 0.02, 0.047 
and 0.08. Remember that the curves on this plot, which correspond to specific values 
of Af2, terminate at the appropriate upper limit on tua (derived from the restriction 
Ml/2 < 800 GeV). This limit decreases as A^j increases. 

Requiring Vt^ h"^ to be confined in the cosmologically allowed range 0.09 — 0.22, one 
finds that A^^ is restricted between 0.02 and 0.08. Note that the upper limit on Afj 
does not depend on our restriction on Mi/2- On the contrary, the lower limit on Afj 
is somewhat dependent on the particular choice one makes for this restriction. This 
deserves further study which would require going beyond the simplifying assumption of 
a common supersymmetry threshold Ms- It should be pointed out that this lower bound 
on Af2 is anyway evaded if there exist additional contributions to the cold dark matter 
of the universe from particle species other than x- 

Fig.^ shows the cosmologically allowed region in the itla — Afj plane obtained from 
the above considerations. We see that tua can vary only between about 95 and 216 GeV. 
The lower (upper) boundary of this region corresponds to Vt^ h? = 0.09 (0.22). The left 
boundary corresponds to M1/2 = 800 GeV (0.09 <Vl^h'^< 0.22). Along this line, is 
essentially constant and acquires its maximal allowed value ~ 354 GeV (see Fig.^). The 
minimal value of the LSP mass is obtained at the lower left corner of this allowed region. 
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where Af,^ ^ 0.047, and is equal to about 222 GeV (see Fig.0). We, thus, see that the LSP 
mass ranges between ^ 222 and 354 GeV. The f2 mass is bounded between about 232 and 
369 GeV which makes f2 a phenomenologically interesting charged sparticle. The upper 
(lower) bound corresponds to the upper right (lower left) corner of the region in Fig.^. 
Actually, the whole sparticle mass spectrum is strongly restricted by our considerations. 
Note, however, that the upper bounds on the sparticle masses depend on our choice for 
the maximal allowed Mi/2- This requires a detailed study with inclusion of one- loop and 
supersymmetry threshold effects which may not be negligible for higher M1/2 's. 

The relative contributions cr{ij)/d'eff {{ij) = (xx); (x^2), (7^2X2*'*)) of the three inclusive 
(co)annihilation reactions to o"e// are given in Fig.^ as functions of rriA for the 'central' 
value of = 0.047. The allowed range of niA is 95 — 211 GeV in this case. 

IV. CONCLUSIONS 

We considered the MSSM with gauge and Yukawa coupling unification employing ra- 
diative electroweak symmetry breaking with universal boundary conditions from gravity- 
mediated supersymmetry breaking. We calculated the relic density of the LSP (an almost 
pure bino). Coannihilation of the LSP with the NLSP (the lightest stau) is crucial for 
reducing its relic density to an acceptable level. Compatibility with the mixed or the 
pure cold (with a nonzero cosmological constant) dark matter scenarios for structure 
formation in the universe requires that the lightest stau mass is about 2 — 8% larger than 
the bino mass. This combined with the fact that the LSP mass is restricted to be greater 
than about 222 GeV allows the lightest stau mass to be as low as 232 GeV. 

ACKNOWLEDGMENTS 

We thank B. Ananthanarayan, C. Boehm, M. Drees, N. Fornengo, S. Khalil, A. B. La- 
hanas, K. Olive and N. D. Vlachos for discussions. One of us (C. P.) thanks P. Porfyriadis 
for computational help and the Greek State Scholarship Institution (I. K. Y) for support. 
This work was supported by the European Union under TMR contract No. ERBFMRX- 
CT96-0090 and the Greek Government research grant PENED/95 K.A.1795. 



18 



REFERENCES 



[1] E. Gawiser and J. Silk, Science 280 (1998) 1405. 

[2] Q. Shafi and F. W. Stecker, Phys. Rev. Lett. 53 (1984) 1292. For a recent review 
and other references see Q. Shafi and R. K. Schaefer, |hep-ph/9612478 . 



[3] T. Kajiata, talk given at the XVIIIth International Conference on Neutrino Physics 
and Astrophysics (Neutrino'98), Takayama, Japan, 4-9 June, 1998. 

[4] G. Lazarides, Phys. Lett. B452 (1999) 227. 

[5] A. E. Evrard, MNRAS 292 (1997) 289. 

[6] S. Buries and D. Tytler, Astrophys. Journal 499 (1998) 699; ibid. 507 (1998) 732. 

[7] B. Schmidt et ai, Astrophys. Journal 507 (1998) 46; A. Riess et ai, Astron. J. 
116 (1998) 1009; S. Perlmutter et ai, Bull. Am. Astron. Soc. 29 (1997) 1351. 

L. Krauss and M. S. Turner, Gen. Rel. Grav. 27 (1995) 1137; J. P. Ostriker and P. 
J. Steinhardt, Nature 377 (1995) 600; A. R. Liddle et ai, MNRAS 282 (1996) 281. 



[9] J. R. Primack and M. A. K. Gross, |astro-ph/98l020l 



[10] G. Lazarides and N. D. Vlachos, |hep-ph/9903511| (to appear in Phys. Lett. B). 



[11] H. Goldberg, Phys. Rev. Lett. 50 (1983) 1419. 

[12] J. Ellis, J. S. Hagelin, D. V. Nanopoulos, K. A. Ohve and M. Srednicki, Nucl. Phys. 
B238 (1984) 453. 



[13] M. Drees and M. M. Nojiri, Phys. Rev. D47 (1993) 376; M. Drees, |hep-ph/9703260 



[14] J. Ellis, T. Falk and K. A. Olive, Phys. Lett. B444 (1998) 367; J. Elhs, T. Falk, 
G. Ganis, K. A. Olive and M. Schmitt, Phys. Rev. D58 (1998) 095002; J. Elhs, T. 
Falk, K. A. Olive and M. Srednicki, |hep-ph/9905igT. 



[15] A. B. Lahanas, D. V. Nanopoulos and V. C. Spanos, |hep-ph/ 9906394 . 



[16] K. Griest and D. Seckel, Phys. Rev. D43 (1991) 3191. 



19 



[17] S. Mizuta and M. Yamaguchi, Phys. Lett. B298 (1993) 120; P. Gondolo and J. 
Edsjo, Phys. Rev. D56 (1997) 1879. 

[18] B. Ananthanarayan, G. Lazarides and Q. Shafi, Phys. Rev. D44 (1991) 1613. For 
a more recent update see U. Sarid, [hep-ph/9610341| . 

[19] B. Ananthanarayan, G. Lazarides and Q. Shafi, Phys. Lett. B300 (1993) 245; L. 
J. Hall, R. Rattazzi and U. Sarid, Phys. Rev. D50 (1994) 7048; V. Barger, M. S. 
Berger and P. Ohmann, Phys. Rev. D47 (1993) 1093; P. Langacker and N. Polonsky, 
Phys. Rev. D49 (1994) 1454. 

[20] M. Carena, M. Olechowski, S. Pokorski and C. E. M. Wagner, Nucl. Phys. 
B426 (1994) 269. 

[21] B. Ananthanarayan, Q. Shafi and X. M. Wang, Phys. Rev. D50 (1994) 5980. 

[22] V. Barger, M. S. Berger and P. Ohmann, Phys. Rev. D49 (1994) 4908; M. Drees 
and M. M. Nojiri, Nucl. Phys. B369 (1992) 54; M. Olechowski and S. Pokorski, 
Nucl. Phys. B404 (1993) 590. 

[23] G. Gamberini, G. Ridolfi and F. Zwirner, Nucl. Phys. B331 (1990) 331. 

[24] D. M. Pierce, J. A. Bagger, K. T. Matchev and R. Zhang, Nucl. Phys. B491 (1997) 3. 

[25] H. Baer, M. Brhlik, D. Castaiio and X. Tata, Phys. Rev. D58 (1998) 015007. 

[26] This problem arises in the simplest realization of the 5*0(10) model with universal 
boundary conditions and Yukawa coupling unification, where the electroweak higgses 
exclusively belong to a 5*0(10) 10-plet. However, in a complete 50(10) model of this 
variety, which correctly incorporates fermion masses and mixings, the electroweak 
higgses will receive smaller contributions from other representations too. Under some 
circumstances (see, e. g., G. Lazarides and Q. Shafi, Nucl. Phys. B350 (1991) 179), 
these contributions can affect only the down quark and lepton masses leaving, in 
particular, the top quark mass unaltered. The details of the scheme can then be 
arranged so that mf,(m^) becomes experimentally acceptable, while m,- and tan/3 
are minimally affected. Alternatively, small GUT threshold corrections to gauge (or 



20 



Yukawa) coupling unification can also help to reduce rrih in specific models without 
abandoning universal boundary conditions. An example with a small negative GUT 
correction to and acceptable value of rrif, can be found in T. Blazek and S. Raby, 
Phys. Rev. D59 (1999) 095002. It is interesting that such a correction to as is anyway 



needed for gauge coupling unification |p4|. Our LSP relic density calculation is not 
affected in any essential way by these modifications of the simplest scheme. 

[27] H. Baer, C. Chen, M. Drees, F. Paige and X. Tata, Phys. Rev. Lett. 79 (1997) 986. 

[28] M. Drees and M. M. Nojiri, Phys. Rev. D45 (1992) 2482. 

[29] J. Elhs, G. Ridolfi and F. Zwirner, Phys. Lett. B262 (1991) 477. 

[30] E. W. Kolb and M. S. Turner, The Early Universe (Addison- Wesley, Redwood City, 
CA, 1990). 

[31] K. Griest, M. Kamionkowski and M. S. Turner, Phys. Rev. D41 (1990) 3565. 
L. Roszkowski, Phys. Rev. D50 (1994) 4842; J. Wells, |hep-ph/ 9404219 . 



[33] H. Haber and G. L. Kane, Phys. Rep. 117 (1985) 75; J. Gunion and H. Haber, Nucl. 
Phys. B272 (1986) 1, (E) ibid. 402 (1993) 567. 



21 



FIGURES 

I I I I I I I I I I I I I I I I I I I I I I I I I 




200 



100 120 140 160 

m, (GeV) 



180 



200 



220 



FIG. 1. The mass parameters mo and M1/2 as functions of niA for Af^ = 0.02 (solid lines) 
and 0.08 (dashed lines). 
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FIG. 2. The relevant part of the sparticle spectrum as a function of mA for Af^ = 0.047. 
The LSP mass, for Afj = 0.02, is also included (dashed line). 
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FIG. 4. The cosmologically allowed region in the tua — Afj plane, where ^lsp lies in 
the range 0.09 - 0.22. We also take niA > 95 GeV and M1/2 < 800 GeV. 
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FIG. 5. The relative contributions ^{xx)l^eff (solid line), ^{xf'2)l^ejj (dashed line) and 
(j^-^-(,)^/(Te// (dot-dashed line) of the three inclusive (co) annihilation reactions to cie// as func- 
tions of ruA for Af„ = 0.047. 



